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We study the transport of ultra cold atoms in a tight optical lattice. For identical fermions the 
system is insulating under an external force while for bosonic atoms it is conducting. This reflects 
the different coUisional properties of the particles and reveals the role of inter-particle collisions in 
establishing a macroscopic transport in a perfectly periodic potential. Also in the case of fermions 
we can induce a transport by creating a coUisional regime through the addition of bosons. We 
investigate the transport as a function of the coUisional rate and we observe a transition from a 
regime in which the mobility increases with increasing coUisional rate to one in which it decreases. 
We compare our data with a theoretical model for electron transport in solids introduced by Esaki 
and Tsu. 

PACS numbers: 05.60.-k, 03.75.Ss, 03.75.Lm 
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The motion of particles in periodic potentials is the un- 
derlying process for fundamental transport phenomena 
like electric current in metals. Quantum mechanically, 
particles in a periodic potential can be described with 
Bloch states. Without an external force, the particles 
can move freely through the potential and in the case of 
a non-interacting sample the system acts like a perfect 
conductor. Under a constant external force, the periodic 
potential is tilted and the new stationary states are lo- 
calized Wannier-Stark states 0. In the absence of inter- 
actions the particles cannot change their quantum state 
and the latter system behaves like an insulator for DC 
currents. Instead, in the presence of interactions colli- 
sions can change the quantum state of the particles and a 
macroscopic current is established. At the onset of inter- 
actions an increasing coUisional rate is therefore expected 
to favor a current through the potential whereas at high 
coUisional rate the current is hindered by collisions. The 
latter regime is well known from solids where scattering 
with phonons and impurities provide an extremely large 
coUisional rate and the conductivity decreases linearly 
with increasing coUisional rates 0- However the limit 
of low coUisonal rate, where the role of collisions is re- 
versed, is experimentally not accessible in solids. With 
the development of semi-conductor superlattices j3l , this 
regime could be entered and phenomena like negative 
electric conductivity could be observed 0, 0, Q) but a 
completely non-interacting system is not achievable even 
in superlattices. 

In this work we use ultra cold atoms in an optical lat- 
tice to investigate the transport in periodic potentials 
induced by an external force starting from the limit of 
zero interaction. Such kind of systems have already been 
used with success to study solid state phenomena like the 
Wannier-Stark ladder or Josephson junctions 0. Here 
we take advantage of the unique possibility of controlling 
both the scattering process and the parameters of the lat- 
tice to study the transition from an ideal insulator to a 
real conductor in a perfectly periodic potential. The use 



of indistinguishable fermionic atoms allows us to create 
a completely non-interacting system since inter-particle 
collisions are forbidden at ultra low temperatures. On 
the other hand we can add a scattering mechanism in 
a controllable way by using a mixture of fermions and 
bosons or pure bosonic samples. 

For the experiments we employ a mixture of ultra cold 
fermionic '^"K and bosonic ®^Rb atoms in a magnetic 
trapping potential. The spin polarized fermions are sym- 
pathetically cooled by forced evaporation of the bosons 
1^ llflj. For potassium (rubidium) atoms, the magnetic 
trap has an axial and radial oscillation frequency of 
= 27r X 24(16) s^i and lu, ^ 2tt x 280(190) s^^. The 
experiments were carried out with samples at tempera- 
tures between 300 and 400 nK. The number fermions can 
be varied between 2 x 10** and 10^ which corresponds to a 
Fermi temperature of 300 — 400 nK. By changing the final 
ramp of the radio frequency evaporation, we can adjust 
the number of bosons in the mixture and with a sweep 
below the trap bottom it is also possible to remove the 
bosons. The temperature of the mixture is always above 
the critical temperature for Bose-Einstein condensation. 
During the last 500 ms of the evaporation ramp, the opti- 
cal lattice is switched on adiabatically allowing for ther- 
malization within the lattice. The lattice is formed by 
two counter-propagating laser beams and is aligned in the 
horizontal plane along the weak axis of the magnetic trap. 
Its potential is given by U{x) = Uo/2{l — cos(47ra;/A)), 
where A = 830 nm is the wavelength of the laser. The 
lattice depth is measured in units of the recoil energy 
Uq = sEr with = h?k'^/2m and k = 2tt/X. In the 
axial direction, the combined periodic and harmonic po- 
tential has two different kinds of eigenstates that belong 
to the first energy band. On both sides of the potential 
localized stationary states develop for sufficiently high 
magnetic field gradients. These states are the analog to 
the Wannier-Stark states in the case of a linear poten- 
tial and typically extend over a few ten lattice sites. In 
the trap center, the states extend symmetrically around 
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FIG. 1: (a) Sketch of the experimental sequence (b) Evolution 
of the center of mass (CM) position of a cloud of fermions. A 
pure fermionic sample (circles) does not move to the trap cen- 
ter, whereas an identical sample with an admixture of bosons 
reveals a current through the potential (squares). The data 
are fitted with a sum of an exponential decay and an ini- 
tial damped oscillation as described in the text (continuous 
lines) . The expansion time of the cloud is 8 ms and the lattice 
height is s = 3. The temperature and the atom number of 
the fermions are T = 300 nK and = 5 X lO"*, the number of 
admixed bosons is Nb = 1 x 10^. 



the potential minimum and resemble harmonic oscillator 
eigenstates which are modulated by the periodic poten- 
tial The temperature of the samples is comparable 
to the recoil energy and is chosen in order to have a sig- 
nificant occupation of the localized states. The typical 
1/e^-radius of the cloud in the direction of the lattice is 
100 //m corresponding to roughly 250 lattice sites. In the 
two radial directions, the atoms occupy the radial har- 
monic oscillator states. To study the transport of the 
particles along the lattice, the magnetic trap is suddenly 
shifted ^3 ill the direction of the lattice by a fraction of 
the extension of the cloud (displacement Xd). The har- 
monic confinement acts like an external driving force and 
we monitor the center of mass (CM) position of the cloud 
(see Fig.CJi). 

The evolution of the CM motion for an initial displace- 
ment of Xd =35 fim is shown in Fig.QjD. The open circles 
show the motion of a pure fermionic sample. The filled 
squares show the evolution of the fermionic sample in the 
presence of bosons. After an initial damped oscillation, 
which is attributed to particles in the harmonic oscilla- 
tor like states, the pure fermionic cloud remains displaced 
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FIG. 2: Decay time r of a cloud of fermions in a mixture 
with bosons in dependence on the collisional rate (dots). The 
first data point at zero collisional rate was taken for a pure 
fermionic sample. The number of fermions is A = 50 000 with 
a temperature of 350 uK. The number of bosons was changed 
from 25 000 to 300 000 corresponding to a change in the inter- 
species collisional rate between 40 s~^ and 550 s~^ p^ . The 
lattice height for the two species was sk = 3 and SRb ~ 9, 
the initial displacement was Xd ~ 35 nm. The solid line is a 
drift time, calculated from Eq. Q for a linear potential with 
Bloch oscillation frequency lubo = 2n x 35 s~^ (see text). 



in the trapping potential. This is due to the asymmet- 
ric occupation of the localized states after the displace- 
ment which gives rise to an offset of the cloud with re- 
spect to the equilibrium position. Without collisions, the 
fermions cannot change the population of the single par- 
ticle states. Because Landau-Zener tunnelling into higher 
bands is also negligible for our parameters, the offset can- 
not vanish and the system is insulating. With an ad- 
mixture of bosons however, the fermions rapidly move 
towards the equilibrium position [l^ . This macroscopic 
transport corresponds to a DC current. To quantify this 
current, we fit an exponential decay to the long time tail 
of the data. For the fermions in the mixture we find 
a decay time of t = 260 ± 30 ms, whereas for the pure 
fermionic sample the decay time is longer than 5 s which 
is comparable to the lifetime of the atoms in the opti- 
cal potential. This experiment proofs that in a perfect 
lattice interactions between the particles are needed to 
establish a macroscopic current under an external force. 

To investigate the role of collisions in more detail, we 
have changed the collisional rate by changing the num- 
ber of bosons in the mixture. Fig.|21shows the decay time 
for the fermionic cloud in a range where the collisional 
rate was changed over more than one order of magni- 
tude. Starting from the coUisionless regime, we observe 
a decrease of the decay time with increasing collisional 
rate. This is what one expects if the collisions assist the 
hopping between different localized states. For high colli- 
sional rates, the experimental data show a slight increase 
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of the decay time with increasing coUisional rate. In this 
regime the number of bosons is much higher than the 
number of fermions and the bosons can be regarded as 
a thermal bath for the fermions. The fermions exhibit a 
drift motion and the coUisions with the bosons impede 
the current through the potential like in a electric con- 
ductor. This change in the mobility of the current carri- 
ers is also known from negative differential conductivity 
in semiconductor superlattices 001 (NDC), which is re- 
lated to the phenomenon observed here. In the regime of 
NDC, the current through the potential decreases when 
the applied voltage is increased. This is due to the tighter 
localization of the electron wave function which reduces 
the transition probability of a hopping event between the 
localized states 0, In our experiment, we do not 
change the transition probability but the rate of transi- 
tion inducing collisions. Despite the different mechanism 
for NDC only the product of the transition probability 
and the coUisional rate determines the final hopping rate. 
This formal identity allows us to compare our experimen- 
tal data with the theoretical model that was introduced 
by Esaki and Tsu |3] to describe NDC. The authors calcu- 
late the drift velocity of electrons in a periodic potential 
under a constant external force. They introduce a phe- 
nomenological scattering rate 7 and show that the drift 
velocity depends on the ratio of the Bloch oscillation fre- 
quency in the linear potential ujbo and the scattering rate 
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with vq = XAE/h being the tunnelling speed through the 
potential and AE being the width of the first band. A di- 
rect adaptation of the above equation to our dynamics is 
rather complicated because we have a spatially varying 
Bloch oscillation frequency [l^ and an inhomogeneous 
system. However, we can compare the initial velocity 
Vi of the center of mass that we observe in the experi- 
ment with the drift velocity calculated from Eq. ^ for a 
uniform system in a linear potential. To determine the 
Bloch oscillation frequency in this potential, we take the 
force that initially acts on the center of mass after the 
displacement and we identify the scattering rate 7 with 
the average coUisional rate between the fermions and the 
bosons. Because the initial velocity of the center of mass 
is connected to the decay time hy t = Xd/vi we can also 
compare the decay time r with the inverse of the drift 
velocity Vd- The result is shown in Fig.|21 where we have 
computed the solid line leaving vq in |^ as a free fit- 
ting parameter. In spite of the simplification the model 
reproduces well both the drop of the decay time at low 
coUisional rates and its slight increase for high coUisional 
rates. This supports the phenomenological interpretation 
given above. 

Because we have to deal with inter-particle collisions, 
our damping mechanism is different from that in super- 
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FIG. 3: Decay time of a cloud of bosons for different lattice 
heiglits. Tlie initial displacement was a::d ~ 10 fim. The con- 
tinuous line is an exponential fit to the data. The exponent 
is given by e~"^^'^ . 



lattices and we have to ask if the assumptions made in 
Ref. to derive expression ^ are fully valid in our sys- 
tem. The phenomenological scattering rate 7 describes 
dissipative scattering processes, where the electrons can 
arbitrarily exchange energy and momentum with an ex- 
ternal thermal bath. In our system, no energy exchange 
with the lattice is possible because the lattice is free of im- 
purities or excitations and momentum can only be trans- 
ferred to the lattice in multiples of the Bragg momentum 
via umklapp scattering processes. As mentioned above, 
in the case of a mixture of fermions and bosons, the lat- 
ter might be regarded as a thermal bath for high atom 
numbers and a dissipative scattering channel exists. For 
small numbers of admixed bosons and for pure bosonic 
samples, where the assumption of having a thermal bath 
is questionable, we find the same phenomenology. This 
indicates that also in this case a dissipative mechanism is 
still present, possibly related to the coupling to the two 
radial degrees of freedom. 

We have also studied how the transport of the atoms is 
affected by different lattice potentials which can be eas- 
ily tuned by changing the laser intensity. In experiments 
carried out with thermal bosons, we have measured the 
dependence of the decay time on the lattice height, such 
as shown in Fig.O The data show a rapid increase of 
the decay time with increasing lattice height. This can 
be explained with a reduction of the tunnelling proba- 
bility between neighboring lattice sites with increasing 
lattice height. In expression (Q) the tunnelling speed ap- 
pears as a scaling factor for the drift velocity. Even if 
one takes into account a spatially varying Bloch oscilla- 
tion frequency ujbo {x) and an inhomogeneous scattering 
rate 7(2;), the role of vq does not change. Thus, we can 
write for each single particle a differential equation of 
the form x/vq = f{ujBo{x),j{x)), whose solution scales 
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in the time domain with vq- Consequently, also the be- 
havior of the center of mass scales with vq and the decay 
time must be proportional to the inverse of the band- 
width. For a sinusoidal potential, the bandwidth can be 
expressed in terms of Mathieu functions and we find that 
for s < 10 the bandwidth is well described (the maximum 
error is smaller than 10 percent) with an exponential drop 
of the form AE = E^e~''^^'^. One therefore expects an 
exponential increase of the decay time with increasing 
lattice height. The exponential fit in Fig.|21 demonstrates 
well that this dependence is accomplished. We find a 
numerical value for the factor in the exponent of the fit 
of 1.6. For other experimental data sets with different 
temperatures and initial displacements we derive values 
ranging from 1.5 to 4.5. 

We can now identify two crucial processes that are 
needed for the macroscopic transport through a periodic 
potential in the presence of an external force. The first 
one is the tunnelling from one lattice site to the next 
one. However, the coherent nature of the tunnelling pro- 
cess leads to a localization of the particle. Therefore an 
additional dissipative process is needed to destroy the lo- 
calization of the particle wave function. If one of these 
two mechanisms is missing, the system is insulating, as 
we observe it for non-interacting fermions and in the limit 
of deep lattices. 

In conclusion we have demonstrated that under an ex- 
ternal force the macroscopic transport of particles in a 
perfect periodic potential requires an interaction between 
the particles. Indeed, in the non-interacting case, we find 
an insulating behavior of the system. In the interact- 
ing regime, we observe a transport, which is mediated 
by collisions. We have investigated the dependence of 
the transport velocity on the coUisional rate and on the 
lattice height. A comparison with a semiclassical band 
model introduced for electrons in superlattices reveals 
a good qualitative agreement although the microscopic 
dissipative mechanism is different. In this first study of 
the transport of fcrmionic atoms in optical lattices the 
interaction was provided by an admixture of bosons. In 
the future it will be interesting to prepare the fermions 
in two different spin states and to extend the study to 
the regime of high quantum degeneracy. An even more 
flexible way of tuning the interaction can be achieved by 
using a Feshbach resonance. Such a system of interact- 
ing fermions would be interesting for possible studies of 
the influence of fermionic superfluidity on the transport 
along a periodic potential. 
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